Landau- Zener Tunnelling in a Nonlinear Three-level System 
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We present a comprehensive analysis of the Landau-Zener tunnelling of a nonlinear three-level 
system in a linearly sweeping external field. We find the presence of nonzero tunnelling probability 
in the adiabatic limit (i.e., very slowly sweeping field) even for the situation that the nonlinear term 
is very small and the energy levels keep the same topological structure as that of linear case. In 
particular, the tunnelling is irregular with showing an unresolved sensitivity on the sweeping rate. 
For the case of fast-sweeping fields, we derive an analytic expression for the tunnelling probability 
with stationary phase approximation and show that the nonlinearity can dramatically influence the 
tunnelling probability when the nonlinear "internal field" resonate with the external field. We also 
discuss the asymmetry of the tunnelling probability induced by the nonlinearity. Physics behind 
the above phenomena is revealed and possible application of our model to triple-well trapped Bose- 
Einstein condensate is discussed. 
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I. INTRODUCTION 

Avoided crossing of energy levels is a universal phe- 
nomenon for the quantum non-integrable systems where 
the symmetry break leads to the splitting of degener- 
ate energy levels forming a tiny energy gap. Around the 
avoided crossing point of the two levels the Landau-Zener 
tunnelling (LZT) model provides an effective description 
for the tunnelling dynamics under assumption that the 
energy bias of two levels undergoes a linear change with 
time|l|]. It is a basic model in quantum mechanics and 
has versatile applications in quantum chemistry Jjj, col- 
lision theory [3| , and more recently in the spin tunnelling 
of nanomagnets 0,0, Bose-Einstein condensates (BEC) 
and quantum computing pj, to name only a few. 

LZT model has been extended to many versions tak- 
ing diverse physical conditions into account: LZT prob- 
lem with a time- varied sweeping rate|8|, LZT model with 
a fast noise from the outer environment , LZT model 
with periodic modulation [T(i El, and so on. Among 
them, LZT in a nonlinear two-level system is one of 
most interestin g m odels and attracts much attention 
recently [l2l IHi. Il4. fl5| . In this model, the level ener- 
gies depend on the occupation of the levels, may arise in 
a meanfield treatment of a many-body system where the 
particles predominantly occupy two energy levels. The 
nonlinear LZT model not only demonstrate many novel 
behavior of great interest in theory but also has impor- 
tant applications in spin tunnelling of nanomagnets [l6l | 
and a Bose-Einstein condensate in a double-well poten- 
tial or in an optical lattice 0,0. However, 
since most of the problems of interests involve more than 
two energy levels, with transitions between several levels 
happening simultaneously 0, Il9l l20t pll . for example, 
BECs trapped in multiple wells [221 123 . Y2M. l2q . spin tun- 
nelling of nanomagnets with large spin, etc. It is natu- 
rally desirable to extend the above nonlinear tunnelling 



to the multi-level situation. 



In present paper, we consider the simplest multi-level 
system — three-level system, to investigate its compli- 
cated tunnelling dynamics in the presence of nonlinear- 
ity. Because quantum transitions may happen between 
several levels simultaneously, the LZT in the nonlinear 
three-level model show many striking properties distin- 
guished from that of the two-level case. In the adiabatic 
limit we will show that, for a very small nonlinear pa- 
rameter that the energy levels still keep the same topo- 
logical structure as its linear counterpart, the adiabatic- 
ity breaks down manifesting the presence of a nonzero 
tunnelling probability. This is quite different from the 
two-level case, where the break down of the adiabaticity 
is certainly accompanied by a topological change on the 
energy levels. More interestingly, the tunnelling is irregu- 
lar with showing an unresolved sensitivity on the sweep- 
ing rate, a phenomenon attributed to the existence of 
chaotic state. In the sudden limit, we derive an analytic 
expression for the tunnelling probability under station- 
ary phase approximation and show that the nonlinear- 
ity can dramatically influence the tunnelling probability 
at the resonance between the nonlinear "internal field" 
and the external field. We also discuss the asymmetry 
of the tunnelling probability induced by the nonlinearity. 
The physical mechanism behind these phenomena is re- 
vealed and possible application of our model to triple- well 
trapped Bose-Einstein condensate is discussed. 



The paper is organized as follows. In Sec. II we intro- 
duce our nonlinear three-level LZT model and calculate 
its adiabatic levels. Section III discusses LZT among 
the levels. Section IV gives a possible application of the 
model to the triple- well trapped BEC. 



2 



II. THE MODEL AND ADIABATIC LEVELS 



We consider following dimensionless Schrodinger equa- 
tion 
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with the Hamiltonian given by 
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where v is the coupling constant between the neighboring 
levels; c is the nonlinear parameter; the energy bias 7 is 
supposed to be adjusted by a linearly external filed, i.e., 
7 = at, a is the sweeping rate; oi, ct2 ; a 3 is probability 
amplitude in each level and the total probability |ai| 2 + 
l a 2| + |c*3 1 is conserved and set to be unit. 

When the nonlinear parameter vanishes, our model re- 
duces to the linear case and the adiabatic energy levels 
£(7) = 0, ±^y/j 2 + 2v 2 (Fig. 1(a)) derived by diagonal- 
izing the Hamiltonian (2). Tunnelling probability T nm 
(n, 7Ti = 1,2,3) is defined as the occupation probability 
on the m-th level at 7 — > +00 for the state initially on 
the n-th level at 7 — > —00. For the linear case, the above 
system is solvable analytically and the tunnelling proba- 
bilities can be explicitly expressed as jU 



1 — exp — — — 
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FIG. 1: Adiabatic three levels at v = 1.0. (a) linear case, (b) 
weak nonlinearity case of c = 0.1. (c) strong nonlinearity case 
of c = 10.0. 



TTV 

T 12 = 2cxp(- — 
2a 



I , TTV 



(4) 



S 2 e 



isxt '."I , a 2 = \/l - sx — S2e Ma * , a 3 = 
troduce the relative phase Ox = ai — 9 a2 , 6> 2 = a3 — 8 a2 . 
In terms of sx, Ox and S2, 62, the nonlinear three- level sys- 
tem is casted into a classical Hamiltonian system, 
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The others are T 21 = T 23 = T 32 = r 12 ,r 31 = r 13 , 
Tn due to the symmetry of the levels. 

With the presence of the nonlinear terms, we want to 
know how the tunnelling dynamics in the above system 
is affected. In our discussions, the coupling parameter is 
set to be unit, i.e., v = 1. Therefore, weak nonlinear case 
and strong nonlinear case mean that c << 1 and c >> 1, 
respectively. As to the external fields, we will consider 
three cases, namely, adiabatic limit, sudden limit, and 
moderate case, corresponding to a << 1, a » 1 and 
a ~ 1, respectively. 

Similar to the linear case, we need to analyze the adi- 
abatic levels of the nonlinear model first. With ax = 



Sx - S 2 )Sx COS 01 



v y(l - sx - s 2 )s 2 cos 9 2 

si,9x and s 2 ,# 2 are two pairs of canonically conjugate 
variables of the classical Hamiltonian system. The fixed 
points of the nonlinear classical Hamiltonian correspond 
to the eigenstates of the nonlinear three-level system, and 
are given by the following equations: 



s\ = -vy/ (1 - sx - s 2 )si sin 0i 
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S2 = -v\/ (1 - si - s 2 )s 2 sin 2 
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By solving the equations (JSJl-JHJ the eigenstates of the 
system are obtained. Accordingly, the eigenenergy is ob- 
tained by e = H e i.e., the energy levels are gained as 
shown in Fig.l. 

For weak nonlinearity, the levels' structure is similar 
to its linear counterpart (fig. 1(b)). For strong nonlin- 
earity (fig. 1(c)), in the mid-level a double-loop topologi- 
cal structure emerges and in the upper-level a butterfly 
structure appears. Because of these topological distor- 
tions on the energy levels, we expect that the tunnelling 
dynamics will dramatically change. 



III. LANDAU-ZENER TUNNELLING 
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In this section we study LZT in the nonlinear three- 
level system both numerically and analytically. First, we 
consider two limit cases: adiabatic limit, sudden limit, 
respectively. Then we will discuss the tunnelling prob- 
ability in general case and investigate the symmetry of 
the tunnelling probability. 



FIG. 2: The tunnelling probability Tn, r 22 , r 33 (full circles) 
as functions of a for different nonlinear parameter at v — 1.0. 
The dash lines represent the results from the linear Landau- 
Zener model for comparison. 



A. Adiabatic limit (a <C 1) 

In adiabatic limit, the characters of the tunnelling 
probabilities should be entirely determined by the topol- 
ogy of the energy levels and the eigenstates' properties 
(corresponding to the stability of the fixed points in clas- 
sical Hamiltonian system), according to the adiabatic 
theorem [2£||22]]. So, we expect that, for the weak nonlin- 
earity case, an initial state started from any levels (upper, 
mid or lower) will follow the levels and evolves adiabat- 
ically, as a result, no quantum transition between levels 
occurs; for the strong nonlinearity, an initial state from 
the lower level is expected to evolve adiabatically keeping 
stay on the ground state, leading to zero adiabatic tun- 
nelling probability, whereas for the state initially from 
the mid or upper level, due to the topological change 
of the level, it can not move smoothly from left-side to 
the right-side. Transition to other levels happens at the 
tip of the loop or butterfly. Consequently, the adiabatic 
tunnelling probability is expected to be nonzero. 

However, the above picture is only partly corroborated 
by our directly solving the Schrodinger equation using 
forth-fifth order Runge-Kutta adaptive-step algorithm, 
as shown in Fig. 2. 

On the one hand, Fig. 2 clearly shows that, for the 
strong nonlinearity case, as we expect, no tunnelling for 



the state from the lower level, but a serious adiabatic tun- 
nelling is observed for the states from upper two levels. 
In particular, we find that the tunnelling probability as a 
function of the sweeping rate shows an irregular oscilla- 
tion. We associate this irregularity to the chaotic state. 
To demonstrate it, we plot in Fig. 3 the Poincare section 
of the trajectories for c = 10 before and after the tip of 
the butterfly structure of the upper level in Fig.lc. It 
shows that, before the tip, the eigenstate corresponds to 
the fixed point surrounded by quasi-periodic orbit, there- 
fore is stable. As the state evolves to the right tip of the 
butterfly, it contact with chaotic sea, after that the state 
become chaotic. The characteristics of the chaos is sensi- 
tive on the parameters, therefore the chaotic state is re- 
sponsible for the irregular tunnelling probability exposed 
by Fig.2h,i. 

On the other hand, Fig. 2 also shows that for the weak 
nonlinearity, even though the adiabatic levels keeps the 
same topological structure as the linear case, there is 
still nonzero tunnelling probability for the state started 
from the mid-level. The tunnelling also shows some kind 
of irregularity. This phenomenon counter to our naive 
conjecture from observing the topological structure of the 
adiabatic levels. 

To explain this unusual phenomenon, we need make 
detailed analysis on the property of the fixed points of 
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FIG. 3: Poincare section of the trajectories for c = 10 before 
and after the tip of the butterfly structure of the upper level 
in Fig.lc. 



the classical system Hamiltonian(7), corresponding to the 
eigenstates of the mid-level. 

We plot quantity S\ as the function of 7 in Fig. 4 (a,b), 
we see the adiabatic evolution of the eigenstate breaks 
down around 7 = —2 due to the nonlinearity (Fig. 4 (b)). 
This adiabaticity breakage is caused by the change on the 
property of the fixed point corresponding to the eigen- 
state of the mid-level. This is revealed by investigating 
the Hamiltonian-Jaccobi matrix obatined by linearizing 
the nonlinear equations (f5f)— (fill) at fixed points, 
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(12) 



We solve the eigenvalues of Hj for different 7 and plot 
our results in Fig. 4c. These eigenvalues can be real, com- 
plex or pure imaginary. Only pure imaginary eigenvalues 
corresponds to the stable fixed point, others indicate the 
unstable ones. In Fig. 4 (c),we can see the eigenvalues 
are complex number (i.e., their real parts are not zero) 
around 7 = 0,±2. The corresponding fixed points are 
unstable. For other regions, the eigenvalues of Hj are 
pure imaginary. Therefore, even though no topological 
structure changes on the level structures, the instability 
of the fixed point corresponding to the mid-level leads to 
the breakdown of the adiabaticity manifesting the irreg- 
ular nonzero tunnelling probability exposed by Fig.2e in 
the adiabatic limit. 

The above instability mechanism occurs for any 
smaller nonlinear perturbation. Let us make some an- 
alytic deduction as follows. Note that the fixed points of 
equations l|5|l- Hll|) can be accurately calculated if c = : 
s l = s 2 = 2+^- 0? = 0, 6»£ = tt for 7 > 0, and (9? = tt, 
#2 = for 7 < 0. By employing the perturbation theory 



FIG. 4: The variety of si with 7 when the eigenstate (0, 1, 0) 
evolve adiabatically at v = 1.0, a = 0.0001. (a) Linear case, 
(b) Nonlinear case at c = 0.1. (c) shows the real parts of the 
eigenvalues of Hj. 



using c as small parameter, we can get the fixed points for 



small c : s? 



(1-7 2 ) 



5TC7, 



'1 — 2+r 1 ~ 4(2+7 a )° °2 — ' 4(2+7^)' 

6± = 0, 82 = tt for nonlinear case. Substituting them 
into equation (|12f) . we can obtain the eigenvalues of Hj 
by solving the following quartic equation: 



(64+12807 4 )a: 4 +(64+c z +13447' i )a;^+(16+c z +3527^) = 

The useful quadratic discriminant is A = 40967 4 — 
2432c 2 7 2 + (c 4 -128c 2 ). I n linear case, c = 0, A = 40967 4 
is always larger than zero, which means that the solutions 
for x are pure imaginary, thus the fixed points are sta- 
ble. For small c, lim A < 0, the real part of the solutions 

~ c/16, while the imaginary part ~ y/2/2. As a result, 
the fixed point corresponding to mid-level becomes un- 
stable around 7 = for any small nonlinearity, implying 
the break down of the adiabatic evolution of states on 
the mid-level. 



B. Sudden limit (a > 1) 

The sudden limit corresponds to nonadiabatic LZT. 
The tunnelling probability does not relate much to the 
structure of the levels. In this limit a weak nonlinear- 
ity does not affect the tunnelling probability, however, a 
strong nonlinearity can dramatically influence the tun- 
nelling dynamics. 

In this limit, we can derive the analytical expression 
of the tunnelling probabilities using the stationary phase 
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approximation (SPA). As a demonstration, wc concen- 
trate on the mid-level, i.e. to calculate 1^2 which equals 
to 1 — T 2 i — r 23 . Because of the large sweeping rate a, a 
quantum state would stay on the initial level most of the 
time. Thus the amplitudes a\ and a 3 in the Schrodinger 
equation (1) remain small and | ct2 1 ~ 1 all the time. 
A perturbation treatment of the problem becomes ad- 
equate. 

We begin with the variable transformation, 



a>i = a 1 exp| — i 



(1 + 1 
\2 4 
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= a' 2 exp[-i J (| |a 2 | 2 ) dt], (14) 



a 3 = a 3 exp[— i 



-l + l\a 3 \ 2 )dt]. (15) 



As a result, the diagonal terms in Hamiltonian arc trans- 
formed away, and the evolution equations of a 1 , a 2 , 



a 3 become: 
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We need to calculate the above integrals self- 
consistcntly. Due to the large a, the nonlinear term in 
the exponent generally gives a rapid phase oscillation, 
which makes the integral small. The dominant contri- 
bution comes from the stationary point to of the phase 
around which we have 

«;=-j /^^/(i+f |ai|2 ~£H- (i6) 
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This expression can be differentiated and evaluated its 
result at time to- A standard Fresnel integral with the 



result 



dt 



-I to 



is obtained. Combining this 
with the relation (15), we come to a closed equation for 



ai 



a 3c 
2 + T 
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For given a, c, v, a\ in the above equation can be ob- 
tained. The tunnelling probability 
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The alliance of Eq.(20)and (21) gives the analytic expres- 
sion on the tunnelling probability T 23 in the sudden limit. 
Compared with our numerical simulation it shows good 
agreement at c/v < 130, c/v > 130 a clear deviation is 
observable (Fig. 5a). It is due to the resonance between 
the "internal field" and the external field leads to the 
invalidity of our assumption |a 2 | ~ 1, as we show latter. 

Similarly, to calculate T 21 , we consider following equa- 
tion, 
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dt exp[i 
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Differently, in this case we may have three stationary 
phase points that are solutions of equation (24) when c < 



27 V tF^u 2- ' but only one solution otherwise, as demon- 
strated in Fig. 6. We denote them as a 3 i, a 3 2, a 33 from 
smallest to largest. For small c, a 3 i is around —a/2 , and 
the other two solutions locates at the two sides of the ori- 
gin. In this case, we simply take a 3 = a 3 \ + a 32 + a 33 . 



Then 



— J* dtexp[ij*^(t-to) 2 dt]. (25) 



The tunnelling probability 

r 2 i 



,2 

a 3| +QO 
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The alliance of the Eq. (24,26) will give the approximate 
solution of the r 2 i. Compared with our numerical simu- 
lation it shows a good agreement at c/v < 110, whereas 
for c/v > 110 a clear deviation is observed (Fig. 5b). 
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FIG. 5: Comparison between our analytic results using SPA 
(full circles and crosses) and the numerical integration of the 
Schordinger equation (1) (solid lines). A cross is used to de- 
noted the invalidation of SPA 
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What happens around c/v = 110 that leads to the 
break down of our stationary phase approximation? The 
reason is the resonance between the "internal field" and 
the external field. Let us recall the exponent in the in- 
tegrand of Eq.(22), we find the effective sweeping rate 
should be the difference between the change rate of the 
"internal field " (i.e., 1 0.3 1 ) and the sweeping rate of the 
external field. At c/v — 110, we find the two frequen- 
cies become almost identical, leading to the invalidity of 
SPA assumption of rapid phase oscillation. This reso- 
nance accompanied by the bifurcation of the stationary 
phase points. Crossing c/v — 110 we observe the num- 
ber stationary phase points changes from three to one, as 
shown in Fig. 6. The resonance breaks the SPA leading 



to serious transition from level 2 to level 1 , consequently, 
at c/v > 130, our assumption |a 2 | <~ 1 become invalid, 
and our approximation on the r 2 3 from SPA is no longer 
good as shown in Fig. 5a. 



a/V 




FIG. 7: The contour plot of tunnelling probability P22 as the 
functions of the scaled sweeping rate and nonlinearity. 



C. General property of the nonlinear tunnelling 
probability 

The nonlinear tunnelling probability as the function 
of the two scaled quantities a/v 2 and c/v, show many 
unusual properties. Taking mid-level tunnelling r 22 for 
example, we make a large numerical exploration for a 
wide range of parameters, to demonstrate the general 
property of the nonlinear tunnelling probability in Fig. 7. 
In general, increasing the sweeping rate will reduce the 
probability of tunnelling to upper or lower level and the 
positive nonlinearity usually suppresses the probability 
of the state's staying in the mid-level, because that the 
nonlinearity with positive c can be regarded as a kind of 
repulsive potential. This repulsive self-interaction make 
particle tend to transition to lower level more easily, and 
this transition becomes more serious at the occurrence 
of the resonance between the "internal field" and exter- 
nal field. The occurrence of the resonance is clearly ex- 
posed by the boundary between the white regime and 
dark regime in Fig. 7. In the white regime, due to the 
resonance, the nonlinearity dramatically changes the tun- 
nelling probability. 

The other issue we want to address is the symme- 
try. The nonlinearity makes levels deform and there- 
fore break the symmetry between upper level and lower 
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FIG. 8: r nm as the function of a for c = 10 (open pentacles), 
c = 40 (open circles) at v — 1.0. Dashed line denotes the 
linear case for comparison. 



level, consequently, the relations r 2 i = r 23 = r 3 2 = 
r 12 ,r 31 — r 13 ,r 33 — Tn hold in linear case break in 
presence of the nonlinearity. For our three-level system, 
the symmetry breaking is clearly exposed by Fig. 8 of 
showing the tunnelling probability T nm as the functions 
of a/v 2 for c — 10, c = 40. In the linear case, we have 
T 2 i = T2 3 = r 3 2 = Ti2, however, with the presence of 
the nonlinear, Fi2,r2i increases whereas the r2 3 ,r 32 de- 
creases. The similar things happen for the r 3 i,ri 3 and 
r33) Til- The above symmetry breaking may be observed 
experimentally 28J . 



IV. CONCLUSION AND APPLICATION 

In conclusion, we have made a comprehensive analysis 
of the Landau-Zcncr tunnelling in a nonlinear three-level 



system, both analytically and numerically. Many novel 
tunnelling properties are demonstrated and behind dy- 
namical mechanism is revealed. 

Our model can be directly applied to the triple-well 
trapped BEC and to explain the tunnelling dynamics 
between the traps [24l l_!5j . In a triple-trap v(r), a 
BEC is described by Gross-Pitaevskii equation (GPE) 

ifiWgQ. = „£v 2 *(r,i)+ [v(r)+g \*(r,t)\ 2 

under the mean-field approximation, where go — 47rh aN , 
m the atomic mass and a the scattering length of the 
atom-atom interaction. The wave function W(r, t) of 
GPE is the superposition of three wave functions de- 
scribing the condensate in each trap^^i i- e -; ^( r >t) = 
ipi(t)4>i(r) + ip2(t)fa(r) + ip3(t)(j) 3 (r). When we study 
the tunnelling of three weakly coupled BEC in traps 1, 
2 and 3, the dynamics of the system is described by the 
nonlinear Schodinger equation with the Hamiltonian, 



H = 



(27) 

where E° a = /(^ |V0 Q | 2 + v (r) \<j) a \ 2 )dr {a = 1,2,3) is 

the ground state energy for each trap. c a — J go |0 Q | 4 dr 
(a = 1,2,3) stands for atom-atom interaction, i.e., non- 
linear parameter. if 12 = — J (y-V^iV^2 + ^(r)^i0 2 )c?r 
is the coupling matrix element between trap 1 and 2. 
K 2 3 = -/(^V^> 2 V(/>3 + v{r)(j)2^)dr is the coupling 
matrix element between trap 2 and 3. For simplicity, we 
only consider the case that these two coupling matrix el- 
ements are the same and there is no coupling between 
trap 1 and 3, i.e., K12 = K23 = K, K13 = 0. The en- 
ergy bias can be adjusted by tilting the trapping well and 
the nonlinearity can be adjusted by the Feshbach reso- 
nance technique. We hope our theory will stimulate the 
experiment in the direction. 
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